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We derive the representation of the nonequilibrium steady-state distribution function 
which is expressed in terms of the excess free energy production. This representation re- 
sembles the one derived recently by Komatsu and Nakagawa [Phys. Rev. Lett. 100 (2008), 
030601] resting on the use of microscopic time-reversal symmetry, but our representation 
applies also to sheared granular systems in which such a symmetry is broken. 

§1. Introduction 

Developing statistical mechanics for nonequilibrium steady states is one of the 
most challenging problems in theoretical physics .^^^^ Among the most remarkable 
outcomes from those studies devoted to such a development have been the generalized 
Green-Kubo relatiorP and various forms of work and fluctuation theorems .'^''^ A 
suggestive representation of the nonequilibrium steady-state distribution function 
in terms of the excess entropy production has recently been developedjIS'l^ and its 
connection to the steady-state thermodynamics has been discussed.'^ 

However, most of these studies rest on the use of microscopic time reversality or 
the local detailed balance, and the outcomes therefrom cannot be applied literally to 
macroscopic, irreversible dissipative systems like granular fluid^^ despite manifest 
similarities.'^ Recently, we have demonstrated that the generalized Green-Kubo 
relation and the integral form of the fluctuation theorem can be derived without 
assuming the microscopic time reversality,'^ thus extending the major achievements 
so far in the nonequilibrium statistical mechanics to irreversible dissipative systems. 
A purpose of the present paper is to further continue this line of work, and to explore 
the representation of the steady-state distribution function which resembles the one 
derived in Refs.Ej) and [9]), but applies also to irreversible dissipative systems. This 
will be exemplified for uniformly sheared granular systems. 

The paper is organized as follows. In SecEJ after introducing basic equations of 
motion for uniformly sheared granular systems, we summarize exact relations which 
serve a basis in our subsequent formulation. We then derive in Sec. [3] the represen- 
tation of the steady-state distribution function which resembles the one developed 
in Refs. |8]) and[9l), but is applicable also to systems in which the microscopic time- 
reversal symmetry is broken. The paper is summarized in Sec. [H Appendix |X] is 
devoted to a technical derivation of the result which is used in the main text. 
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§2. Exact starting points 

In this section, we introduce basic equations of motion for sheared granular 
systems and corresponding Liouville equations. We then summarize some exact 
relations which serve a basis in our subsequent formulation. 

2.1. SLLOD equations 

The system of our interest consists of N smooth granular particles of mass m 
which is immersed in a volume V and is subjected to stationary shearing char- 
acterized by the shear-rate tensor k. We assume that each granular particle is a 
soft-shere, and the contact force acts only on the normal direction. For a simple 
uniform shear with velocity along the x-axis and its gradient along the y-axis, which 
we consider throughout this paper, the shear-rate tensor is Kq,^ = 'jdax^iSy with 7 
denoting the shear rate. Under homogeneous shear, a linear streaming velocity pro- 
file ft • r is induced at position r. Newtonian equations of motion describing such a 
homogeneously sheared system are the SLLOD equations^ 

fi = ^ + K-ri, (2-la) 

m 

p, = Ff+Ft'^-K.p,. (2.1b) 

Here refers to the position of the ith particle, and is the peculiar, or thermal, 
momentum defined with respect to the streaming velocity k ■ ri. The conservative 
force P'f^'^ = —dU/dvi with the total interaction potential U is given by a sum 
F^^'^ = X^jyj F^^j"^ of the elastic repulsive forces exerted on the ith particle by others 

= = &{d - n,)f{d - r,,)f,„ (2-2) 

where u{r) is the pair potential; rij = Vi — Vj, rij = \rij\, Vij = rij/vij] d denotes the 
particle diameter; 0{x) is the Heaviside step function. Typical forms of the elastic 
repulsive force are /(x) oc x (linear model) and f{x) oc x^/-^ (Hertzian model). Sim- 
ilarly, the viscous dissipative force Fj^'^-* due to inelastic collisions between particles 
is represented by a sum F^^^^^ = "Ej^i -f^ij'**^ of two-body contact forces 

Fjf ) = -0{d - r,,)C(d - r,,)ig,^ ■ h,)hr (2-3) 

It is proportional to the relative velocity g^j = fj — ifj = (pj — p,j)/m + k • Vij of 
colliding particles, and for this reason, the SLLOD equations (|2-ip are not invariant 
under the time-reversal map {'"i,Pi,7} — {fi,—Pi-,—"i}- The amount of energy 
dissipation upon inelastic collisions is characterized by the viscous function (^(x) 
which is typically assumed to be constant or modeled as C(x) cx x^/^. 

2.2. The Liouville equation 

For nonequilibrium systems, the form of the Liouville equation commonly used 
for Hamiltonian systems should be properly generalized to take into account the 
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effect of phase-space compression.^' The Liouville equation for phase variables, say 
A{r), where F = {r^ ,p^) stands for a phase-space point, is given by 



dt 



A{r) = r--^A{r)^iC{r)A{r). 



(2-4) 



The operator iC{r) is referred to as the p-Liouvihean. The explicit expression 
for iC{r) corresponding to the SLLOD equations (|2-ip shall be written for later 
convenience in the form 

iC{r) = iCo{r) + iAC{r) (2-5a) 



with 



i4/;(r) = ^[(K.ri).^-(Kp. 
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(2-5b) 
(2-5c) 



Here, iCo{r) represents an "unperturbed" part which is derivable from the Hamil- 
tonian describing the internal energy of the system 



(2.6) 



while the rest iAC{r) is due to the driving forces towards the nonequilihrium states. 
Here, not only the shearing force, but also the viscous dissipative force should be 
considered as the driving force. This notion is clear for the sheared granular system 
under study: even when the shearing force is absent, the system gets of out equilib- 
rium due to the dissipative force describing inelastic collisions between particles. 

The formal solution to the Liouville equation (|2-4p can be written in terms of 
the p-propagator ex.'p[iC{r)t] as 



A{r{t)) = exp[i/:(r)t] A{r). 



(2-7) 



Hereafter, the absence of the argument t implies that associated quantities are eval- 
uated at t = 0, and the dependence on F shall often be dropped for brevity like 
A{t) = A{r{t)) and iC = iC{r). 

On the other hand, the Liouville equation for the phase-space distribution func- 
tion f{r,t), normalized such that f dF f{r,t) = 1, is given bj® 



df{r,t) 

dt 



dr 



rf{r,t) 



r-^+A{r)]f{r,t) 



ic\r)f{r,t). (2- 



The operator iC'^ is called the /-Liouvillean, and A(r) = {d/dF) ■ F is referred to 
as the phase-space compression factor. For the SLLOD equations ()2Tp . one has 



d 



A{r) = y^.F 



-I (vis) 



— ^^Oid-r,,)ad-r, 



(2-9) 
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It follows from ([24]) and ([2^ that 

iC\r) = iC{r) + A{r). (2-io) 

Thus, the Liouville operator is non-Hermitian in the presence of the viscous dissipa- 
tive force. Noticing from ()2-9p that A{r) is completely determined by the viscous 
dissipative force, iC^ shall also be decomposed in the form of (j2-5p as 

ic\r) = iclir) + iAc\r) (2-iia) 

with 

iclir) = iCo{r), iAc\r) = iAC{r) + A{r). (2-iib) 

The formal solution to the Liouville equation ()2-8p reads 

/(r,t) = exp(-i/:tt)/(r,o), (2-12) 



where exp(— i>C^t) is called the /-propagator. For later convenience, the initial dis- 
tribution function f{r,0) shall be chosen as the equilibrium one generated by iCj^, 
i.e., as the solution of 

i4(r)/(r,o) = o. (2-13) 

It follows from iC^ = iCo and ()2-5bp that /(P, 0) for the present system is given by 
the canonical distribution function of the inverse temperature /3, 

/(r,0) = ^f-'""^"^^ with Zil3) = J dre-^"^^\ (2-14) 

As demonstrated in Ref. 13), nonequilibrium steady-state properties based on the 
canonical initial distribution function do not depend on the choice of f3. 

Let us summarize here for later use some relations between /- and p-Liouvilleans 
and corresponding propagators. When iC^{r) acts on the product of the initial 
distribution function f(r,0) and a phase variable A{r), there holds 

z£t(r)[/(r, o)A{r)] = A . [tfA] = a^ ■ [tf] + fr-^A 

= Air)[ic\r)fir,o)] + f{r,o)[iC{r)A{r)]. (2-15) 

One can show that iC and iC^ are adjoint operators: 

j dr [i£{r)A{r)] B{r) = - j drA{r) [iC^r)B{r)]. (2-i6) 

This property can be proved from the integration by parts. By a repeated use of 
this property, the following relation for the propagators can be derived: 

J dr [e^^(-^)*A(r)] B{r) = J drA{r) [e-''^^(-^)*s(r)]. (2-17) 



Representation of the nonequilihrium steady- state distribution function 5 



2.3. Generalized Green-Kubo relation 
With the identity 

e-^'"'^''^' = 1 + f dse-'^'^'^^l-iC^r)], (2-18) 

the formal solution (j2-12p can be expressed as 

/(r, t) = fir, 0) + f ds e-^^'^^> [-i£t (r)/(r, o)] 

Jo 

^/(r,o)+ rd5e-^'(^)^[/(r,o)^2(r)], (2-19) 

Jo 

where in the final equality we have defined the dissipation function n{r). From an 
explicit calculation, one finds for the present systercPJ 

n{r) = -fi^cj^y{r) - 2pn{r) - A{r). (2-20a) 

Here, the shear-stress tensor aafs{r) and Rayleigh's dissipation function TZ{r) are 
respectively given by 

a^,{r) = + rUFif + f/J))] , (2-20b) 

^(^) = i E E - ''^^■)^('^ - ■ h,f. (2-20c) 

i j^i 

Let us summarize here some important properties concerning the dissipation 
function Q{r). First, one can show that the equilibrium average of the dissipation 
function is zeroff^) 

(i7(r))=0. (2-21) 

Here and in the following, the ensemble average (• • • ) is defined over the initial 
equilibrium distribution function /(f ,0), i.e., 

(•••)^ j drfir,0)--- . (2-22) 

Second, there hold^^ 

dr{t)f{r{t),o) = dr f{r,o)e-^o^'^^^^'^\ (2-23) 

This relation is a consequence of the conservation of the number of ensemble members 
within a comoving phase volume.^^ Finally, we notice that, if we define the time- 
dependent free energy via f3F{t) = j3{H{t)) — S{t) in terms of the Gibbs entropy 
S{t) = -/dr/(r,t) log/(r,t), one finds pT{t) = {Q{t)) since H = -^a^y - 27^ 
and S{t) = {A{t)). Thus, the average dissipation function divided by /3 can be 
identified as the free energy production rate. 
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The nonequilibrium ensemble average {A{t)) shall be defined via 

(Ait)) = J drf{r,Q)A{t) = j drf{r,t)A{Q). (2-24) 

The two representations in terms of f{r,0) or f{r,t) are equivalent because of the 
relation (pTT]) . Substituting (p49]) for f {r,t), one obtains 

{A{t)) = (A(0)) + r ds {A{s)m), (2-25) 



Jo 

in deriving which we have used ()2T7p . 

The system is said to be in a nonequilibrium steady state if the ensemble averages 
of all phase variables become time- independent. Let us notice that the long-time 
limit of (j2-25p becomes constant for systems that exhibit mixingpl' which is to be 
assumed in the following. This feature can be demonstrated as follows. By taking 
the time derivative of (|2-25p . one finds 

= {A{t)m)- (2-26) 

For systems that exhibit mixing, all the long-time correlations between phase 
variables vanish, and we obtain for t — )• oo 

±{A{t))^{A{t)){m)=0, (2-27) 

where we have used the property (|2-2ip . This indicates that the long-time steady 
state average of an arbitrary phase variable becomes constant, i.e., 

lim (Ait)) = (A)ss, (2-28) 

t— >oo 

where the steady-state average, denoted by (• • • )ss hereafter, is obtained from the 
t oo limit of (f2^ : 

/•CXI 

(A),, = (A(0)) + / ds {Ais)niO)). (2-29) 
Jo 

This is called the generalized Green-Kubo relation relating the steady-state average 
to the integral of the transient time-correlation function,SJ'li31i and applies to systems 
arbitrarily far from equilibrium. It reduces to the conventional Green-Kubo relation 
if the shearing force is weak and the viscous dissipative force is neglected. 

§3. Cumulant-expansion representation of the distribution function 

In this section, we derive a representation of the phase-space distribution func- 
tion which resembles the one developed in Refs.H]) and^), but without assuming the 
microscopic time reversality. Our derivation is certainly related, but not identical, 
to the one adopted in Refs.EI andEJ. 
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3.1. Conditioned "averages'' and cumulant- expansion representation 

Let us first notice that, using the relation (|2-17p with r{t) = e*^*i^, the distri- 
bution function fi'jjt) at a phase-space point 7 can be expressed as 

/(7, t) = Jdr fir, t) 6{r -^) = Jdr f{r, o) <^(r(t) - 7). (s-i) 

With this notion, we shah introduce the conditioned "average" {C)^-t in which the 
phase-space point at time t is constrained to 7: 

^ Jdrcf{r,o)6{r{t)--f) _ J drc f{r,o)5{r{t) --f) 

/dr/(r,o)5(r(i)-7) " fh,t) ' ^''"^ 

Notice that {C)-y-t so defined is not really an averaged quantity (this is why we put 
quotation marks) since, our equations of motion being deterministic, it is completely 
determined by the single path in the phase-space that ends up at the point 7 at time 
t. We shall also introduce the conditioned average {C)~f-Q, in which the phase-space 
point at time t = is constrained to 7, by setting t = in ()3-2p : 

^ !drcf{r,o)6{r--f) _ f drc f{r,o)6{r --f) 

jdrfir,o)6ir-^) " 7(7,0) • ^'''^ 

By inserting C = e^-^o d^s^^i-s) (|3.3p _ ^e have 

/(7,o) (ei/orf^«(-«))^.o = J dr f {r, 0)6^0 <i'^(n-s)) s{r-'r) 

= J dr f{r, 0) ^(-^("^ 6{r - 7). (34) 

By setting r = r{t), one obtains, since r{s) = r{t + s) 

/(7,0)(e^/cN-^(-))^.o = J dr(t)/(r(t),0)el/orf^^(^W)5(r(t)-7). (3-5) 
Using (1223]), this leads to 

/(7,0) (e^i'o'>'««(~^))^.o = J (ir/(r,0)e-^i'o'^^«('^W) 5(r(t)-7) 

= /(7,t)(e-^i'o"^^(^))^;t, (3-6) 



where the second equality follows from the definition ()3-2p . We therefore obtain the 
following symmetrical representation of the distribution function: 



/(7,t) (eUo'^^^(-))^. 



7;0 



(3-7) 



In the following, we shall manipulate the right-hand side of ()3-7p using the 
cumulant expansion 

00 ^ 

log(e^)^;, = 5]-(y%. (3-8) 

k=l 
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Here (y'^)^.^ denotes the fcth-order cumulant of Y defined by 



iY%;t^Q^log{exp[uY])^,, ^. (3-9) 



u=0 

For example, {Y)'^.^ = {Y)j-t and {Y'^)'^.^ = {Y'^).y-t — {Y)j-f It is useful to note that 
for any constant yo there holds 

{{Y-yo)%.^t = {Y%.^, for k>2, (3-10) 



which can be derived on the basis of (j3-9p . 

By applying the cumulant expansion to ()3-7p . we obtain 



where we have introduced 

0-= [ dsf2{-s), 9+= [ dsf2{s). (3-12) 
Jo Jo 

Notice that the t dependence is dropped from 0± for notational simplicity. 

We would like to express ()3-lip in terms of excess quantities. For this purpose, 
one needs to introduce averages, so let us define (only formally for the moment) 

/}„ = lim - [ (f?(-s))^.o, ^2+= lim - [ ds {n{s))~..t. (3-13) 

t— >oo t Jq t^oo t Jq 

Under what circumstances these limits exist (not in the mathematical sense but in 
the physical sense) will be discussed below. We will also argue below that these 
averages are independent of the choice of 7, and this is why we dropped the possible 
dependence on 7 from the notation D±. Keeping these facets in mind, we shall 
introduce the excess quantities 

01^ = / ds[f?(-s) - /}_], 0;^= [ ds[n{s)-D+]. (3-14) 
Jo Jo 

In view of the notice below ()2-23p . 0^? divided by /3 can be identified as the excess 
free energy production. With the application of (|3-lUp . there hold 

{6%.^, = {{e'^J')%,o, (0+)^;* = {{0T)%t ioT k>2 (3-15) 

for those cumulants in the second term on the right-hand side of ()3-lip . For the first- 
order cumulants appearing there, we simply exploit the relation 0± = 0'^ + tD±. 
In this way, we obtain from ()3-lip 

00 ^ 

+ E ^ {((0i^)')7;o - i-^)H{0T)%;t} , (3-16) 

fe=2 

which is the starting point of the following formulation. 
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3.2. Expansion in the small "degree of nonequilihrium" 

The representation (|3-16p is formal, and will be useful in practice only if it is 
dominated by a first few terms in the expansion. We shall therefore explore in the 
following the form of the distribution function by restricting ourselves to the small 
"degree of nonequilihrium" to be characterized by a small parameter e. We then 
examine the order in e of each term on the right-hand side of (|3-16p . The usefulness 
of such an expansion has been demonstrated in Ref. [TO]) . where the distribution 
function valid to the second order in e has been exploited to derive thermodynamic 
relations for nonequilihrium steady states. 

As we have noticed above, not only the shearing force, but also the viscous 
dissipative force should be considered as the driving force towards nonequilihrium 
states. The magnitude of the shearing force is determined by the shear rate 7, 
whereas from ()2-3p that of the viscous dissipative force is characterized by Q/m 
assuming a constant viscous function, C{x) = Q. The small degree of nonequilihrium 
shall therefore be characterized by small 7 and C/™-; ^'^d the parameter e shall be 
chosen as the larger of 7 and C/m so that the correction iAC to the unperturbed 
Liouvillean i£o in (|2-5p reads iAC = 0(e), i.e., 

iC = iCQ + iAC with iAC = 0{e). (3-17) 

Notice that our characterization of the small degree of nonequilihrium makes sense 
since small C,/m (i.e., quasielastic limit) is involved there. At the same time, we have 
to assume small 7 since our primary interest here is in the transient dynamics occur- 
ring before the steady state is reached. These points are discussed in Appendix [Bl 
We also note that e so chosen has the dimension of inverse of time. This is because, as 
we will see below, e appears mostly in the form of the product er with the relaxation 
time r. It is therefore more preferable to have er dimensionless. 
Similarly to (fSTTp . we have from ([M]) and (|2TT]1 

iC'^ = iCl + iAC'^ with iAC'^ = 0(e). (3-18) 

Correspondingly, p- and /-propagators can be expanded as 

e^^* = e*^«* + 0(et), e"^^^* = e'^^o* + 0(et). (3-19) 

It is also clear from (I23D and (fT^ that 

n = 0(e), (3-20) 

which plays an important role in the following. For later convenience, we shall intro- 
duce the notation FQit) = e^'~°*r evolving under the "unperturbed" p-propagator 
e*^o*. Because of (1349]) . there hold 

r(t) = ro(0 + 0(et), A{r{t)) = A{rQ{t)) + 0{A)-0{et). (3-21) 

Since the initial distribution function is chosen as the solution of iCQ{r)f{r, 0) = 
0, it follows from the relation (j2-15p specialized to the Liouvilleans iCQ and iCo that 

icl{r)[f{r,o)A{r)] = f{r,o)[iCo(,r)A{r)]. (3-22) 



10 



S.-H. Chong, M. Otsuki and H. Hayakawa 



By repeated applications of this result, one obtains 

e<int [f(^r^ o)A{r)] = f{r, o) [e*^«(-^)*A(r)], (3-23) 

which is to be used in the following. 

3.3. Order estimate 

To estimate the order in e of each term on the right-hand side of ()3T6p . one 
needs to estimate the order of 0^. For this purpose, we need to know more about 
the averages f2± defined in (|3-13p . Let us start from a related quantity 

lim - [ ds {f2{s))^.o = lim - [ dsf2{j{s)). (3-24) 

t^oo t Jq t->oo t Jq 

The system is expected to settle to a unique nonequilibrium steady state, irrespective 
of the initial phase-space point 7, for times longer than the relaxation time which 
we denote as r. Here, r shall be chosen so that the system reaches the steady state 
for t > T. Then, by writing 

dsf2{-f{s))= / dsn{j{s))+ / dsf2{-f{s)), (3-25) 







the second term on the right-hand side can be estimated as (i — r)(l7)ss with the 
steady-state average (^?)ss- As a result, there holds irrespective of the choice of 7 



1 



lim - / ds (^?(s))^;o = (^2)ss. (3-26) 

Now, we apply this argument to = limt^aoi^ / 1) Jq ds {f2{s))-y-f The differ- 
ence from (|3-24p is that here what is fixed is the phase-space point 7 at time t. To 
take this into account, we introduce the time t' required to reach the phase-space 
point 7 during the steady-state fluctuations. Then, the effect of fixing the end point 
7 at time t will show up only in the time regime t — r'^s <t. Therefore, by writing 

f ds{Q{s))^^t= r ds{Q{s))^,t+ f ^ ds{Q{s))^^t+ f ds{Q{s))^,t, (3-27) 

Jo Jo Jt Jt-r' 

and noticing that the second term on the right-hand side can be estimated as {t — 
t' — r)(l7)ss, there holds irrespective of the choice of 7 

f2+ = lim - / ds {f2{s))^.t = {n)ss = lim {n{t)). (3-28) 

t— >oo t Jq ' i— >oo 

From this argument, it is also clear that the integrand of O*^ = ds [i^{s) — 
is non-negligible only for < s < r and t — t' < s < t. As will be discussed in the 
next subsection, however, we will eventually set the upper limit t of the integral in 
6*!j^ to r. In this case, since = 0(e), there holds 

0?^ = 0(er). (3-29) 
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One can apply a similar argument to Q- and ©^f', assuming that a steady-state 
is reached in the negative time direction (see below concerning this point), to obtain 

Q = lim - / ds{Q{-s))~..Q= lim(I?(-t)), (3-30) 

t— >oo t Jq ' t— >oo 

irrespective of the choice of 7, and 

0^^ = 0(er). (3-31) 

From (|3-29p and (|3-3ip . one understands that the cumulants of A; > 3 appearing 
in the right-hand side of ()3-16p are of order O(e^r^). In addition, one finds that 
the second-order cumulant term {{0'^)'^)^.q — ((0^)^)^.^ is also O(e^T^), and this is 
demonstrated in Appendix |X1 We therefore obtain from ()3-16p 

log = I {(01^)^;o + {0T)tA + ^ + ^+) + Oie'r'). (3-32) 

Now, what is left is the order estimate of the second term, which we write in view 
of (|3^ and (fCT]) as 

D_ + D+= lim + (3-33) 

For systems possessing the time-reversal symmetry, one can show that j7_ -|- j7+ = 0. 
However, this is not the case for systems in which such a symmetry is broken. 

3.4. On the relation between and 17_ 

Here, we shall deal with Dj- for systems in which microscopic time-reversal 
symmetry is broken. The average i7+ = limt_^oo(^(i)) in the positive time direction 
is well defined since the system does reach the steady state for t > t. However, in 
general, the average i7_ = lim.t^oo{^i—t)) defined in the negative time direction 
does exist for the time-irreversible system. This is clear for the sheared granular 
system under study: in the negative time direction, particles "attain" energy upon 
collisions, and the system continuously heats up so that the steady state is never 
reached. However, by considering the small degree e of nonequilihrium, there might 
be a possibility that such a heat up of the system up to the time scale — r in the 
negative direction is still negligible. We shall therefore explore in the following the 
relation between J7_|_ and i7_ that holds for small e and up to the time scale itr. 
This means that, instead of (j3-33p . we will consider 

+ {{n{-t)) + {f2{t))} (3-34) 

t=T 

evaluated at time r just after the steady state is reached. This is the reason why we 
eventually set the upper limit t of the integral in to r, which is assumed above 
in connection with ()3-29p . 

We start from the generalized Green-Kubo relation (j2-25p for (i7(— t)) and 
{f2{t)). Since (I?(0)) = 0, we have 

{f2(-t)) = f ds{n{s)f2{0)) = - [ ds{f2{-s)n{0)), (3-35) 
Jo Jo 
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{f2{t)) = [ ds{Q{s)f2{0)). (3-36) 
Jo 

Using the decomposition ()3-19p for e*^* and noticing [2 = 0(e), one finds for {Q{t)): 

{f2{t)) = f^ds j drf{r,o) [e'^'f2{r)] n{r) 

= ds J drf{r,o) [{e'^°' + o{es)} n{r)] Q{r) 

= ds J drf{r,o) [e'^'^'nir)] Q{r) + o{eh'^). (3-37) 

Similarly, one obtains for {f}[—t)): 

{Q{-t)) = - j^ds j drf{r,o) [e-''^'[2{r)] f2{r) 

= - j^ds j dr f{r, 0) [e-^^«" f2{r)] n{r) + 0{eH^). (3-38) 
Applying first (|2-17p and then (|3-23p to this expression, one finds 

{[2{-t)) = - f ds j drf{r,o)[2{r) [e^^o^i?(r)] +o{eh^) 

= -{n{t)) + 0{eH^). (3-39) 
Thus, from the identification ()3-34p . we have for small e 

t + =0{e^T^), (3-40) 

t=T 

for the time- irreversible system at time r just after the steady state is reached. 

§4. Summary 

Collecting our results so far, we have for the steady-state distribution function 
fssir) for the small degree e of nonequilibrium 

Ur) = fir, 0) exp [ i { (0-)r;o + {0T)r;t} ] + ©(eV^). (4-1) 

As argued above, fss{r) for the system in which time-reversal symmetry is broken 
should be considered as the one just after the system entered the steady state at 
time r, fssir) = fir,t)\t=T- Up to the presence of r in the correction terms, one 
understands that the representation ()4-ip for the steady-state distribution function 
is essentially the same as the one derived in Refs. [8]) and iSj) in that it is expressed 
in terms of the excess free energy productions and that it is valid up to the second 
order in e, i.e., it holds beyond the linear-response regime. 

However, there is a subtle point connected with the presence of r in the cor- 
rection terms. In fact, there holds er = 0(1) on general grounds. This can be 
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understood from the generalized Green-Kubo relation {fi{t)) = Jq ds (f?(s)f?(0)) by 
noticing that the left-hand side is 0(e) whereas the right-hand side is O(e^r) since 
(i7(s)J7(0)) decays to zero on the time scale of r. Thus, all the correction terms in 
the representation ()4-ip are of order 0(1), and hence, it is not a converging repre- 
sentation. In this sense, the nonequilihrium steady-state distribution function might 
not exist. This point will be discussed in more detail in our subsequent publication. 
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Appendix A 

Order estimate of the second-order cumulant term 



In this appendex, we show that 

((01^)')^;o - {{Or)X;t = 0{e'r'). (A-1) 
To this end, we notice that with the remark given above ()3-29p there holds 0'^ = 
0{eT) in leading order for small e. Therefore, what (|A-ip claims is that such a 
leading-order contribution in the left-hand side cancels out. 

Using the notation in ()3-2ip . the leading-order term is given by 

6)^-0= f ds[Q{ro{±s)) - fi±o]- (A-2) 
Jo 

Here, fi±Q is defined as in (j3-13p but with e='=*^'^ replaced by the unperturbed e^*^''^. 
Then, to show the validity of (|A-ip . it suffices to demonstrate 

{{0'-o)X;0 - ((0?o)')^;r = 0. (A-3) 
Using the property (|3-10p and expanding cumulants in terms of moments, this means 
that what we have to show are the equalities 

ds{Q{-s))^.fl= f ds{n{s))^,t, (A-4) 





t ft ft nt 

ds I ds' {Q{-s)Q{-s'))~f.o= ds ds' {f2{s)n{s'))-y.t, (A-5) 
Jo ' Jo Jo 

in which p- and /-Liouvilleans are given by iCo and iC^, respectively. 

We start from deriving ()A-4p . From the definition ()3-3p of the conditioned aver- 
age, the left-hand side of (|A-4p is given by 

{ni-s))^,o = 1^ ds I dr [e-^«(^)^f2(r)] /(r, o)5(r - 7). (a-6) 
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Applying first ()2-17p and then ()3-23p to this expression, we obtain 



ds {f2{-s))^.o 



1^ ds I dr f}{r) f{r, o) 6{ro{s) - 7). (a-7) 



On the other hand, from the definition (|3-2p of the conditioned average, the right- 
hand side of (|A-4p is given by 



ds{n{s)) 



/(7,0) 70 



ds / dr 



g^£o{r).^(P)l f^r, 0)6{ro{t) - 7). (A-8) 



Here we have noticed that /(7,t) = /(7, 0) when the /-Liouvillean is given by iC^. 
Again, by applying first ()2-17p and then p-23p to this expression, one obtains 

ds {n{s))^,t = 1^ ds I dr Q{r)f{r, o)6{ro{t -s)-j) 

= 1^ ds I dr Q{r) f{r, o) 6{ro{s) - 7), (a-9) 

and thus the desired equality ()A-4p is derived. 

One can derive the equality ()A-5p in a similar manner, so we only show the 
course of the derivation. First, the left-hand side of (|A-5p can be manipulated as 



t r-t 

ds I ds' {Q{-s)Q{-s'))^.o 
Jo 

1 

7(7^ 



t rt 

ds / ds' I dr 



JQ 
t rt 



-iCoir)i 



Q{r) 



-iCo{r)s' 



Q{r) f{r,mr-i) 



^—j^dsj^ds'j dr Q{r)Q{ro{s - s'))f{r, 0) 6{r^{s) - 7). (A-io) 



/(7 



For the right-hand side of ()A-5p . one can proceed as follows: 



t rt 

ds / ds' {f2{s)n{s'))^,t 



Jo 

1 

7(7,0) 70 ^"Jo 



t i-t 

ds / ds' I dr 



e'^''^^>n{r) e'^''^^>' n{r) /(r,o)(5(ro(t) -7) 
/(70) Jo I ^(^)^(^o{s' - s))f{r, o)S{ro{t - - 7) 



"7(70)/ "^'L '^'''/ '^^^(^)^(^o(^'-* + ^))^(^'0)'^(^o(s)-7) 

= 1^ ds 1^ ds' I dr f2{r)f2{ro{s - s'))f{r, o)S{ro{s) - 7). (A-ii) 

In this way, the desired equality (|A-5p is derived. 
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Appendix B 

Balance equation for the kinetic temperature 

In this appendix, we shall derive a balance equation for the kinetic temperature 
T{t) = 2{K{t)) / {DN) defined for a D-dimensional system in terms of the average of 
the kinetic energy K = X]jP^/(2m) at time t, and discuss its implication. We first 
notice that the time derivative of K is given by 

^ = Y.--P^ = Y.-- i^f' + ^S^'O -^-T. V^vdm^ (B-1) 

i i i 

where we have used ()2-lbp and introduced the notation A : B = 'Ylia p ^apBpa- 
Using Pi/m = Vi — k ■ Vi from (j2-lap . one obtains 



i i 

= T.-^-Ff' + \Y.O^rFf'-n:cr. (B.2) 



2 

1;3 



In the second equality, we have used g^j = Vi — rj, f\^^^^ = FiY^\ Newton's 

third law -f^™^ ~ ~FiY^\ and the definition (|2-20bp of the shear-stress tensor. 
The first term on the right-hand side can be related to the time derivative U of 
the total potential energy since U = r-j ■ {dU/dvi) and pf'^^ = —dU/dvi. The 
second term can be expressed in terms of Rayleigh's dissipation function since TZ = 
-(1/4) ^ijQij ■ fS/"-" (see ([231) and (12^2031) ). We therefore obtain from ([B^ 

k = -U- ^a^y - 27^, (B-3) 



where we have used the explicit form = ^Sax^py of the shear-rate tensor. 

To simplify the following discussion, we assume a constant viscous function, 
C{x) = C- Then, C, can be factored out from IZ (see ()2-20cp ). and we define ^{t) via 
CT^it) = {TZ{t))/V. Let us also introduce the nonlinear viscosity r]{t) via 7ry(t) = 
— {axy{t)) /V . One then obtains from (|B-3p the following time-evolution equation, or 
the balance equation, for the kinetic temperature T{t): 

rw = -^M*)) + A,v*)-^W). (B.4) 

Here n = N/V denotes the average number density. 

As discussed in connection with ()2-27p , the time derivative of the average of any 
phase-space variable vanishes in the steady state. Thus, there holds from ()B-4p 

Tit) ^ -^i\s - T^C^ss = for t ^ oo, (B-5) 
Vn Un 

since limt^oo(t^ (*)) = 0, where r/ss = limt^oo ??(i) and = limf_j.oo T^it). Thus, the 
steady-state kinetic temperature Tss = Wmt-^oo T{t) is determined by the balance 
between the viscous heating {2j'^r]ss/{Dn)) and the collisional cooling {Ad^'Dgs/iDn)). 
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Let us consider an implication of (|B-5p . As noted below ()2-14p . the steady-state 
averages such as ??ss and are independent of the choice of the initial inverse 
temperature /3, and depend only on the "thermodynamic" parameters (A^, 7) and 
on the system-specific parameters characterizing elastic repulsive force /(x) in ()2-2p 
and the viscous function in (j2-3p . the latter being assumed to be constant 
Cix) = C here. In the following discussion, N, V, and the parameters specifying 
f{x) play no significant role, and the dependence of the steady-state averages on 
them shall be suppressed. We shall also introduce dimensionless quantities, to be 
distinguished with tilde such as 7 and (, which are necessary to properly characterize 
the degree of nonequilibrium of the steady state. This can be done in terms of d, m, 
and T(t) as in Ref. Iisp . or in terms of d, m, and the spring constant k assuming a 
linear spring model for f{x). In any case, one obtains from ()B-5P 

%s(7,C) 

which is an analogue of the relation derived in Ref. |15l) based on the fluctuating 
hydrodynamics. The solution to this equation gives 7(C) that depends only on (. 
This means that the dimensionless shear rate 7 cannot be made small by solely 
controlling the shear rate 7, and to realize the small degree of nonequilibrium which 
is characterized by small 7 and C, one has to consider the quasielastic limit C l.'^ 
In the main text, we are interested not only in the steady state which is char- 
acterized only by C in the above sense, but also in the transient dynamics occurring 
before the steady state is reached. In the transient regime, the above argument does 
not hold, and the dynamics there depends on the shear rate 7 as well. Thus, one 
has to take into account both 7 and C in characterizing the transient dynamics. 
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